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ABSTRACT 


Stoneflies (Eustheniidae, Gripopterygidae) were light-trapped at 24 sites on small to large streams over a 300 
m elevation gradient. Catches were small but 8 species occurred; jackknife estimates of total richness (8-11 
spp.) suggest the species list is reasonably complete. Emergence patterns were weakly seasonal with apparent 
peaks in the warm, wet season. Some species ranged extensively along the elevation/stream size gradient but 
Canonical Correspondence Analysis supports the existence of a compositional gradient from lower, larger 
woodland streams to higher, smaller streams in rainforest. We suggest improvements in survey design for 


investigating stoneflies in tropical streams. 
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INTRODUCTION 

Relative to their temperate zone counterparts, 
tropical Plecoptera are incompletely understood. 
Stoneflies can be numerically (Bispo & Oliveira 2007) 
and functionally (Cheshire et al. 2005) important in 
tropical streams. However, stoneflies appear to be a 
significant exception to the general pattern of higher 
biodiversity in the tropics (Hillebrand 2004). Family 
and generic diversity decrease toward the equator 
(Zwick 2000; Vinson & Hawkins 2003; Palma & 
Figueroa 2008; Fochetti & Tierno de Figueroa 2008). 
Regional and local species diversities are less certain. 
Some tropical genera e.g. Anacroneuria are species- 
rich (Stark 2001) and few local surveys are available. 
Here we describe the stonefly fauna of a catchment in 
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the Wet Tropics Area of northeastern Australia. 
Stream research in the Wet Tropics region 
(Pearson 1994) began as faunistic (Pearson et al. 1986) 
and process investigations at single localities and has 
expanded to include comparative, geographically 
extensive analyses (McKie et al. 2005; Pearson 2005; 
Connolly et al. 2008). Our intermediate scale project 
sampled multiple sites along gradients of elevation 
and stream size within a single catchment. Much 
stream research e.g. Cheshire et al. (2005) emphasizes 
functional roles of immature insects typically 
identified to family or genus. In contrast, we 
collected adults which can be identified reliably to 
the species level. In a potentially diverse and 
incompletely known tropical fauna, this was 
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important. The tradeoff, of course, is number of 
specimens; less numerous adults are less likely to be 
collected at all sites where larvae occur. We describe 
occurrence patterns over space and time and 
statistically evaluate the effect of limited collection 
size on interpretation. This problem may be 
especially acute in tropical streams so we suggest 
tactics to improve surveys such as ours. 


MATERIAL AND METHODS 

Running River originates as Birthday Creek and 
flows westward to the Burdekin River from the 
coastal ranges of northern Queensland near the 
village of Paluma (19° 00’ S, 146° 12’ E), 70 km north 
of Townsville. The catchment lies in Mount Spec 
State Forest. Human impacts include limited logging, 
past tin mining, a few roads and a water supply 
reservoir on a tributary. Highest point in the 
drainage is 1020 m. Bedrock is granite and soils are 
shallow and nutrient poor (Congdon and Herbohn 
1993). Rainforest at higher elevations grades 
downward to sclerophyll woodland at our lowest site 
(640 m). Small rainforest streams are densely shaded, 
becoming more open with increasing size and entry 
into woodland. Stream profiles are stepped with 
lower gradient reaches alternating with steeper 
sections and waterfalls. Water at higher elevations is 
weakly acidic (pH 6.2-6.5) and soft (conductivity 
32uS/cm) (Pearson & Connolly 2000). Climate is wet 
but seasonal (Frith & Frith 1985) with a cool, dry 
winter (June-August) and a warm, wet summer 
(December-March). Stream flow reflects the seasonal 
rainfall but varies considerably between years. Our 
study began late in a three year drought but included 
a cyclonic rain January 31-February 1 1994 (948 mm 
at Paluma weather station). This storm scoured the 
streams but much less severely than did a 1991 
cyclone (R. G. Pearson, pers. comm.). Stream 
temperatures 1993-1994 in Birthday Creek (785 m) 
ranged from monthly minima of 9-11°C (July, 
August) to maxima of 21°C (Nov.-Feb.) (Z. Rosser, 
pers. comm.). 

We collected adult insects using battery-powered 
Rothamsted light traps with 12V fluorescent bulbs 
(Williams 1948; Benson & Pearson 1988) at 24 sites. 
The project was done within the confines of an 
academic leave; collections began 15 October 1993 
and terminated 8 July 1994. The missing months are 
times of low insect activity (Frith & Frith 1985; 


Illiesia — http://www2.pms-]j.si/illiesia/ 


Benson & Pearson 1988) but our description of 
seasonal patterns is incomplete. We began with 
regular sampling of two small rainforest streams 
(Birthday Creek, 820 m; un-named tributary of Camp 
Creek, 760 m) as part of life history studies. By late 
November an array of sites, 640-940 m elevation on 
small to large streams, was in place (n=24). The two 
primary sites were sampled during each 2-4 d 
collecting bout (n=22). The remaining sites were 
sampled 1-8 times. Difficult access to some sites, 
equipment failures and others’ need for the 
equipment kept the design from being balanced but 
the 130 successful collections (overnight operation, 
substantial insect catch) on 45 nights were fairly 
evenly distributed over time (12-18 collections per 
month), elevation and stream size. Larvae were 
collected with a kickscreen at the two primary sites 
October-April and at most other sites in October; 
these incomplete data are briefly mentioned in this 
paper. 

Each site was described by two variables, 
elevation and stream size, derived from 1:50,000 
topographic maps and general notes were made on 
channel form, riparian vegetation and canopy 
closure. Stream size was indexed by link magnitude 
(Shreve 1966) which is the total number of first order 
tributaries to a stream segment. The downstream site 
on Running River is order 5, magnitude 285 and 12 
stream km from the source of Birthday Creek. 
Sheldon (1985) used elevation and magnitude to 
summarize stonefly distributions in a North 
American drainage. Some correlation of elevation 
and magnitude is inevitable since the largest stream 
in a drainage network occurs at the lowest elevation, 
but the correlation (r = -0.75) is especially strong here 
because small, low elevation tributaries were dry 
most of the time and were not sampled. However, all 
mapped streams were counted in link magnitudes. 

Taxonomic nomenclature follows Michaelis and 
Yule (1988) and McLellan (1996). Collections are 
deposited at James Cook University (Townsville) and 
the Australian National Insect Collection (Canberra). 

Results are presented in descriptive graphics. 
However, given the variable collecting effort and the 
small number of stoneflies collected, we applied 
statistical procedures to evaluate the strength and 
generality of our results. For example, how 
successful were we in achieving a reasonably 
complete species list? To answer this question, 
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following guidelines in Brose et al. (2003), we used 
the first order jackknife estimator and its variance 
(Heltshe & Forrester 1983): 


Sest= Sobs + (1 - 1V/C) Si 


where Sest is the estimated total number of species, 
Sobs the number of species actually collected, C the 
number of collection units, and S: the number of 
singleton species each collected in a single unit. The 
jackknife is a bias correction; an intuitive explanation 
for increasing Sobs by Siis that each species collected 
just once likely had an uncollected counterpart. 

We used multiple regressions, both linear and 
polynomial, to evaluate patterns of species richness 
over time and the gradient of elevation and stream 
size (log magnitude). Collection number, as log C, by 
month or site was included as a covariate in all 
regressions. For correlation and regression analyses 
we followed Sokal & Rohlf (1981), Zar (1984) and 
routines in SYSTAT 10.2 (SYSTAT 2002). 

Diversity (species richness) is but one measure of 


community structure. To assess compositional 
patterns over time and space we used Canonical 
Correspondence Analysis (CCA) executed in PC- 
ORD 4 (McCune & Mefford 1999). CCA is a 
multivariate constrained ordination technique which 
identifies the maximal correlation(s) between the 
matrix of species occurrences or abundances and the 
corresponding matrix of environmental variables. 
The regression model incorporated in CCA 
accommodates unimodal responses to environmental 
factors. For the temporal analysis we used 
frequencies of species occurrences ((trap 
nights/species/month)/(total trap nights/month)) with 
time (month) and collection number (log C) in the 
biological and environmental matrices respectively. 
Occurrence (presence vs. absence) at sites and 
elevation, stream size (log magnitude) and 
collections/site (log C) were used in the spatial 
analysis. Species occurring just once were omitted 
from CCA and June (no stoneflies) and July were 
combined in the temporal CCA. Statistical 
significance was set at P = 0.05 for all analyses. 


Table 1. Composition and abundance of stoneflies in light trap catches from Running River. 





Individuals Collections Sites 
Eustheniidae 
Cosmioperla wongoonoo tropica (Theischinger) 8 6 2 
Gripopterygidae 
Dinotoperla cardaleae Theischinger 291 41 17 
Illiesoperla barbara Theischinger 18 10 8 
Illiesoperla cerberus Theischinger 4 4 4 
Illiesoperla franzeni (Perkins) 14 4 2 
Kirrama abolos Theischinger 1 1 1 
Leptoperla alata Theischinger 4 4 2 
Leptoperla commoni Theischinger 5 3 3 





RESULTS 

Stoneflies were neither abundant nor diverse in 
light trap catches from the Running River drainage 
(Table 1). Of the 130 collections, only 55 (42%) 
contained stoneflies from 19 sites (79%). The eight 
species collected were dominated by Dinotoperla 
cardaleae constituting 84% of the 345 stoneflies 
collected and found at 89% of occupied sites. In 
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contrast, a single Kirrama abolos was collected. The 
remaining species, although few in numbers, were 
captured in multiple collections at multiple sites. 
Applying the jackknife estimator to the occurrence 
by collection (trap night) (Table 1) yields Sest = 8.98 
with 95% confidence limits (8.00, 10.94) indicating 
that few or no additional species remained 
undetected. 
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Fig. 1. Seasonal occurrence by months of Running River stoneflies in light traps. Frequency = (traps 


occupied/total traps per month) x 100. 


Flight periods were prolonged (Fig. 1). Adults, 
including those of the abundant D. cardaleae, were 
collected October-July. However, some seasonality 
was suggested. The onset of cool, dry weather in 
April reduced captures and no stoneflies were 
captured in June. Stonefly diversity peaked 
January-March. Single species also showed some 
evidence of seasonality. Cosmioperla wongoonoo 
emerged toward the end of the dry season and 


larvae with darkened wing pads were not taken 
after January although Theischinger (1983) 
recorded adults, from various locations, most 
months of the year. The frequency of occurrence of 
D. cardaleae peaked February-March and numbers 
increased also; the six largest collections all 
occurred in March. The high flows generated by 
Cyclone Sadie Jan.31-Feb. 1 did not suppress 
subsequent stonefly captures. 


Table 2. Summary statistics for regressions of species richness (S) on time (t), elevation, 
stream magnitude and number of collections (C). Entries under variables are standard 
partial regression coefficients with probabilities in (). 





Dependent t t? 
a) log (S+1) -0.093 -0.419 
P (0.900) (0.600) 
elevation log magnitude 
b) log (S+1) 0.358 0.314 
P (0.190) (0.205) 


log C adj R? 
0.820 0.706 
(0.001) (0.015) 
log C adj R? 
0.868 0.527 
(0.001) (0.001) 
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Fig. 2. Occurrence of Running River stoneflies at 24 sites over the elevation gradient. Species names are 
sequenced from those occurring at low elevation to those found in higher streams. 
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Fig. 3. Occurrences of Running River stoneflies over the gradient of stream size where magnitude 1 streams 
are the smallest. Species names are sequenced from those occurring in small streams to those found in large 
streams. Overlapping points have been displaced slightly for clarity. 


We present spatial patterns relative to both 
elevation and stream magnitude (Fig. 2, 3). Some 
species (D. cardaleae, Illiesoperla barbara) were broadly 
distributed along both gradients. One species 
(Illiesoperla franzeni) occurred in larger, lower streams 
but four species (C. wongoonoo, Leptoperla spp., K. 
abolos ) characterized smaller, higher streams. Note 
that Leptoperla commoni, occupying 200 m of the 
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elevation gradient, was confined to small streams. 
Conversely, the three Illiesoperla spp. appeared more 
closely associated with elevation than stream size. 
Records were sparse for most species and additional 
collecting undoubtedly will extend and clarify 
habitat useage. For example, in October 1993 we 
made a special effort to locate populations of C. 
wongoonoo but found larvae only at the two sites 
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Table 3. Summary statistics for canonical correspondence analysis (CCA) of a) species 
composition (frequency of occurrence relative to time and number of collections, and b) 
species composition (presence-absence) relative to elevation, stream magnitude and number 
of collections. Coefficients for species and environmental variables for (b) are shown in Fig. 





4. 

a) Time Axis 1 Axis 2 
Eigenvalue 0.187 0.128 
% variance 18.0 12.4 
P 0.508 0.056 
Spp-Envir correlation 0.688 0.807 
P 0.848 0.124 

b) Space Axis 1 Axis 2 Axis 3 
Eigenvalue 0.434 0.291 0.099 
% variance 23.2 15.6 5.3 
P 0.006 0.002 0.012 
Spp-Envir correlation 0.851 0.899 0.554 
P 0.046 0.002 0.084 





where adults were collected and in a larger stream 
(magnitude 22, elevation 785 m), thus supporting the 
elevation preference but extending the range of 
stream size for this species. 

We tested the diversity pattern suggested in Fig. 1 
by regressing log (S+1) on t, t? and log C. The full 
polynomial model (Table 2) was significant (P= 0.015, 
adjusted R=0.71) but only log C contributed 
significantly. Thus the apparent diversity peak in 
January to March was not supported. The CCA of 
species composition (Fig 1) leads to a similar 
conclusion (Table 3). Both t and log C contributed to 
both canonical axes but only the second axis 
approached significance and neither of the species- 
environment correlations was significant. Seasonal 
succession suggested by Fig. 1 was not supported. 
The power of these analyses is low since data were 
grouped into just ten monthly periods. 

To examine spatial trends in richness (0-5 
spp./site) over all 24 sites we regressed log (S+1) on 
elevation, log magnitude and log C. The full model 
(adjusted R?= 0.53) was significant (P=0.001) (Table 2) 
but only log C contributed significantly (P=0.001). 
Thus there were no detectable trends in species 
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richness over the gradient of elevation and stream 
size. Only collecting effort mattered. Note that this 
result was strongly influenced by the five sites 
without stoneflies and few collections (median C=1). 

Unlike richness, species composition varied along 
the elevation-stream size gradient although number 
of collections remained a confounding factor (Table 
3, Fig.4). The primary gradient in Fig. 4 is elevation- 
stream size and the contribution of log C comes 
primarily from the occurrence, probably fortuitous, 
of C. wongoonoo at the two primary sites selected for 
intensive sampling. All three canonical axes were 
significant and axes 1 and 2 represented significant 
species-environment correlations. Thus the broad 
patterns described in Fig. 2 and Fig. 3 were 
supported. 


DISCUSSION 

The stoneflies of Running River included only 
two of four families of Australian § stoneflies; 
Austroperlidae and Notonemouridae are restricted to 
temperate Australia (Michaelis & Yule 1988). The 
observed richness of eight species seems low to 
workers familiar with more diverse regional and 
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Fig. 4. Biplot of Canonical Correspondence Analysis (CCA). Points are species with abbreviations derived 
from names in Table 1. Lines show direction and relative importance (length) of environmental factors and 
collecting effort per site in the two-dimensional space jointly defined by species occurrences and gradients of 


elevation, stream size and collecting effort. 


local faunas. However, Australia has only about 190 
species of which 29 are known from north 
Queensland (Michaelis & Yule 1988; Theischinger 
1988, 1991, 1993). Our jackknife extrapolation, with 
95% confidence interval of 8-11 species, suggests that 
few species were missed. Similarly, Pearson et al. 
(1986) reported eight stonefly species from a well- 
studied single site on Yuccabine Creek 100 km north 
from our location and Monteith & Davies (1990) also 
found eight species at multiple sites along a transect 
in the Bellenden Kerr Range. In contrast, in 
temperate Australia Hynes & Hynes (1975) obtained 
12-20 species at nine sites thoroughly sampled for 
larvae plus adults and accumulated 30 species across 
all sites. Available information is in accord with low 
tropical diversity of Plecoptera (Zwick 2000; Palma & 
Figueroa 2008; Fochetti & Tierno de Figueroa 2008) 
and differs from the conclusion of Lake et al. (1994) 
that the total invertebrate fauna of Birthday Creek 
was much more diverse than that of streams in 
temperate Australia. 

Flight periods (Fig. 1) were extended and, at most, 
weakly seasonal. Similar lengthy emergence of 
stoneflies occurs in Brazil (Froehlich 1991) and weak 
seasonality persists near the equator (Zwick 1976). 
Coffman & de la Rosa (1998) compared emergence of 
chironomid Diptera in Costa Rica and Pennsylvania, 
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U.S.A.; emergence periods were 50-60% longer in the 
tropical assemblages. Emergence of Australian 
species tends to be protracted (Hynes & Hynes 1975; 
see Yule 1985 and Lake et al. 1985 for exceptions and 
discussion) and emergence of Running River 
stoneflies was fairly similar to congeners at 37° S. 

Our study has significant limitations, especially 
the lack of collections from mid-July to October. 
Stoneflies are a cool-adapted group (Hynes 1976; 
Michaelis & Yule 1988) and additional species may 
have emerged in winter as they do in temperate 
regions. However, emergence would have to be 
unusually synchronous not to extend into the 
sampled months. The lack of collections from 
temporary streams at lower elevations is important 
because stoneflies tolerate such conditions (Yule 
1985). 

Light traps may not attract all the species present. 
However, adult aquatic insects in the tropics are 
mostly crepuscular or nocturnal (Edmunds & 
Edmunds 1980) and should be vulnerable to light 
traps. Ernst & Stewart (1985) light-trapped 11 of 12 
stoneflies present in summer at a site in central 
U.S.A. and one of us has trapped many (sometimes 
100s) of North American stoneflies difficult to collect 
by other means. However, independent methods 
(beating, searching, emergence trapping) could be 
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used concurrently to validate light trapping, 
especially in the cooler months. If flight activity is 
temperature dependent (Briers et al. 2003), then small 
cool season catches, which also occur in caddisflies 
(Benson & Pearson 1988) and terrestrial insects (Frith 
& Frith 1985) in this region, may be misleading. 
However, Frith & Frith (1985) attributed seasonality 
of insect catches to plant phenology as the proximate 
factor. Benson & Pearson (1988) considered seasonal 
patterns to be driven by larval life histories rather 
than adult behavior and vulnerability. Life history of 
a caddis species in Birthday Creek (Nolen & Pearson 
1992) supports this view; cool season larvae were 
mostly early instars but late instar larvae, pupae and 
emergent adults were concentrated in the warmer 
months. 

Spatial patterns relative to elevation and stream 
size (Fig. 2, 3) were consistent with broader 
zoogeographic patterns. Six species are part of a 
tropical montane fauna extending northward from 
the vicinity of Paluma through northeastern 
Queensland. Some of these species were confined to 
higher, smaller streams but others were more 
generalized. Two species are distributed primarily in 
southeastern Australia. I. franzeni seems to be part of 
an older eastern Australian fauna of larger, lower 
streams which matches its observed distribution in 
Running River. (This is the first record of I. franzeni 
north of the 350 km Paluma-Eungella gap, a 
significant faunal discontinuity for several taxa of 
aquatic insects (Watson & Theischinger 1984; 
Theischinger 2001)). C. wongoonoo is also part of this 
southeastern fauna but montane and, in Running 
River, was associated with tropical northern species 
in higher rainforest streams. These rainforest 
stoneflies may be vulnerable if a warming climate 
displaces habitat and populations upward into the 
limited and fragmented highlands of the Wet Tropics 
Region (Williams et al. 2003; Wilson et al. 2007). 

Statistical support for our conclusions was 
diminished as a consequence of our opportunistic 
design and our decision to sample more sites but less 
intensively. Prior to our work, stoneflies had been 
collected only from a few easily accessible sites in the 
drainage. We elected to explore a wider range of 
habitats but at considerable cost in replication per 
site. Improved designs, based on the results of this 
survey, are suggested below. 

A more general problem is that we captured few 
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adult stoneflies. Overall, our capture rate was 2.65 
stoneflies/trap night. Similarly Zwick (1976) reported 
988 Perlidae from a 9 m? emergence trap checked 
daily for 16 months in an African stream; i.e. 2.05 
stoneflies/trap night. In addition, a single species 
dominated our collections (84%) and Zwick’s (1976) 
emergence data (91%); the remaining species were 
uncommon. Recent experience of one author (ALS) in 
Bornean streams using beating sheet, light traps and 
kick net larval sampling supports the view that 
tropical stoneflies are scarce. 

Our initial objective, taxonomic resolution, was 
attainable only by collecting adults which are the last 
few survivors of a larval cohort. Asynchronous 
emergence will further reduce daily catches, perhaps 
to the point of non-detection at a specific place or 
time. However, protracted emergence will increase 
the odds that a species will be detected somewhere at 
some time which seems to be the case in our data 
since the species list approached completeness. The 
consequences of rarity are variable data with many 
random absences and strong confounding of 
ecological effects with effort (number of collections) 
as implied by all our statistical analyses. The only 
solution is increased effort applied approximately 
equally at all sites and times. 

A possible design for future work in Running 
River or other Wet Tropics streams would use fewer 
sites, perhaps 8-10, distributed over the gradients of 
stream size and elevation. Increased frequency of 
sampling (2-4 times per month) would increase total 
catch, refine temporal coverage and serve as 
insurance against periods of cyclonic rain and other 
unpredictable events. An appealing possibility is to 
use multiple traps (3-5), spaced far enough apart to 
avoid gear interference, at each site and time. This 
procedure also would increase total catch and allow 
explicit calculation of detection probabilities. This 
level of effort appears extreme. However, 
extrapolating from our catch rates and species 
proportions, a design of 8 sites, 4 traps/site and 
weekly sampling would yield 4410 stoneflies of 
which only 705 would be the seven non-Dinotoperla 
species with catches ranging 13-233 individuals 
(median=59). This level of effort might also detect one 
or more of the additional species predicted by the 
jackknife estimate. Alternative collecting methods 
(beating sheet, emergence traps, Malaise traps) 
should be used to determine bias and efficiency of 
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light trapping. With this level of investment, and 
made possible by our species list and other data, 
larval sampling for life-history and population 
studies would be warranted although the number of 
sites might have to be reduced. 

In conclusion, across multiple sites in a single 
tropical catchment, adult stoneflies were numerically 
rare, of low diversity, weakly seasonal and showed 
moderate habitat specialization for elevation and 
stream size. Tropical stonefly assemblages require 
substantial sampling effort to obtain ecological data 
on any but the dominant species. 
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